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Abstract. We study the problem of a Brownian particle diffusing in finite dimensions 
in a potential given by -0 = 0^/2 where (f) is Gaussian random field. Exact results 
for the diffusion constant in the high temperature phase are given in one and two 
dimensions and it is shown to vanish in a power-law fashion at the dynamical 
transition temperature. Our results are confronted with numerical simulations where 
the Gaussian field is constructed, in a standard way, as a sum over random Fourier 
modes. We show that when the number of Fourier modes is finite the low temperature 
diffusion constant becomes non-zero and has an Arrhenius form. Thus we have a simple 
model with a fully understood finite size scaling theory for the dynamical transition. 
In addition we analyse the nature of the anomalous diffusion in the low temperature 
regime and show that the anomalous exponent agrees with that predicted by a trap 
model. 
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1. Introduction 

Systems with quenched disorder are often taken as paradigms for systems exhibiting 
a structural glass transition. The basic physical idea is that for a sufficiently complex 
and frustrated system without quenched disorder, a single particle sees an effectively 
quenched and random potential due to the other particles. At a mean field level there 
exist models where, to all intents and purposes, this analogy becomes exact. One 
can have two models one with quenched disorder and the other without, but highly 
frustrated, which exhibit the same thermodynamics in the high temperature phase and 
the same glass transition at low temperatures jTI. Even if the frustrated non-random 
system possesses a crystalline ground state not shared by the disordered system, this fact 
is practically irrelevant as this state is dynamically never attained. An analogy between 
mean field spin glass models with one step replica symmetry breaking and structural 
glasses has been put forward j21El; and this theory of the glass transition has become 
known as the random first order scenario. Further extensions of the analogy between 
the dynamics of these mean field spin glasses and finite dimensional glass formers have 
provided much insight into glassy dynamics and related phenomena such as ageing and 
effective temperatures One of the simplest models, which has been extensively 
studied, is the so called toy model. Here one takes a particle in a D dimensional space 
with Hamiltonian 



where ■?/' is a quenched random potential taken from some statistical ensemble. The most 
convenient choice is to take to be Gaussian of zero mean with correlation function 



so the field is statistically isotropic and invariant by translation in space. The statics of 
this model can be solved in the Gaussian variational approximation in finite dimensions 
This approximation becomes exact in limit of the number of spatial dimensions D 
going to infinity if A scales as 



In this limit when the correlation function of F is short range a structural glass transition 
is found in Ej. First, at a temperature there is a dynamical transition, which 
is independent of the precise nature of the dynamics as long as it is local. This 
dynamical transition appears due to the appearance of an exponentially large number 
of metastable states and it can be located via static arguments, thus explaining its 
independence of the precise form of the local dynamics [S]. This transition is not 
accompanied by a strict thermodynamic transition. The thermodynamic transition 
appears at a lower temperature Tg where the entropy of the system becomes very small. 
In the free case, where the harmonic confining potential is removed (the case /x = 0), 
a dynamical transiting occurs at a non-zero temperature, but the static transition is 
occurs only at T = 0. The existence of the dynamical transition in the free case 




(1) 



(V^(x)^(y)) = A(|x-y|) 



(2) 




(3) 
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shows up as a transition between a stationary dynamics characterised by the time 
translational invariance of the correlation function and a form (for unbounded systems) 
of the fluctuation dissipation theorem [Sj and a low temperature ageing regime exhibiting 
ageing in the correlation function and modified fluctuation dissipation relations 0- In 
addition at high temperature the late time behaviour of a Brownian tracer particle in 
the force field generated by the potential is normal and is characterised by a non-zero 
diffusion constant which vanishes at [H]. A little thought convinces one that this is 
not possible in finite dimensions and that there is thus something pathological about the 
large D limit as it is taken in this model. When the correlation length of the random field 
is finite (the case where the correlation function of the field is infinite is another matter 
and anomalous diffusion is clearly possible in this case) one can view the system on a 
coarse grained scale at the order of a few tens of this correlation length. We denote this 
length scale by /*. The process can now be viewed as a discrete random walk between 
neighbouring regions with an exponentially distributed time to jump to a neighbouring 
region whose average is given by the Arrhenius law as r = roexp(^). Here AE is the 
energy barrier associated with moving from one region to another. When the average 
value of r, r is finite the system will look like a random walk in the coarse grained 
picture and we find that 



showing that we should have a finite diffusion constant. When ip is Gaussian we expect 
that the energy barriers AE are also Gaussian and thus obtain a value of r which 
behaves as 



which is referred to as a super-Arrhenius behaviour. The dynamical transition is 
however clearly relegated to zero temperature. This argument is backed up by numerical 
simulations, perturbative and renormalisation group calculations [HI EI- In a finite 
system we will always have finite energy barriers and they will always be overcome by 
activated barrier hopping, all be it after very long but finite times. Mean field models 
have diverging energy barriers and it is this divergence which leads to the dynamical 
transition. Indeed the formula Eq. ^ should ultimately become Arrhenius-like for 
systems where the energy barrier are bounded and we will discuss this point in the next 
section on exact results. 

A commonly used paradigm for glassy systems is the trap model jH] where the 
phase space is considered to be made up of a set of traps denoted by i each of depth 
AEi. The simplest version is that where the traps are on a tree-like geometry, as is the 
case in the random energy model [S]. The generalised random energy model has traps 
within traps and is inspired by the Parisi solution for mean field spin glasses where 
the states (corresponding to the bottom of traps) are organised in with an ultra-metric 
structure Pj. The time spent in a trap is exponentially distributed with mean time 
Tj = To exp(/5A£^j) and the distribution of the AE^ induces the distribution of the Tj. 



r 



(4) 




(5) 
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When the disorder averaged value of Tj = r diverges a dynamical transition occurs. 
Forms of the trap model where the traps are located on a finite dimensional lattice have 
been extensively studied [TIH ITT] . Also a non-random microscopic realisation, based on 
number partitioning combinatorics, of the trap model has been found Intuitively 
we expect the trap model picture to be applicable to the problem of a particle diffusing 
in a short range correlated random potential at sufficiently large time and length scales. 
In this paper we show that for the model we study this is indeed the case. 

Going back to the problem of diffusion in a random potential, we consider a model 
where ip = 0^/2 where is Gaussian. Clearly from the arguments above we do expect 
to see a dynamical/glass transition. We assume that the potential energy barriers also 
behave as the energy itself, we thus assume that the statistics of the barrier heights 
behaves as q;0^/2, where a is some positive constant. Therefore average time spent in 
a trap will behave as 

r~roexp(^— ) (6) 

If we take a Gaussian field (p with correlation function 

(0(x)0(y)) = A(x-y) (7) 

with A(0) = 1, then r diverges for f3a > 1, giving a dynamical transition temperature 
Tg = a. Note that in jTHj the problem of a dipole diffusing in a Gaussian electrostatic 
field was considered. Here the effective potential felt by the dipole is ip = — |(V0)^. 
The above argument indicates that a dynamical transition, indicated by a vanishing of 
the diffusion constant, should occur in this model. However in |^ no clear evidence 
for a phase transition was found in numerical simulations in three dimensions. However 
standard methods for generating the Gaussian field used here and in ^3] use a finite 
number of Fourier modes. Indeed we will show that, in the case studied here where 
exact results are possible, there is a finite size scaling in the number of modes A^. This 
finite size scaling smears out the dynamical transition, just as is the case in the standard 
theory of equilibrium second order phase transitions. 

The Langevin of dynamics we shall study of a particle in a potential ip is given by, 

Xt = Vit) - W(Xi) (8) 

where ri{t) is Gaussian white noise with correlation function 

{r],{t)r]j{t')) = 26i,6{t-t'). (9) 

This choice of white noise amplitude ensures that the stationary measure is the Gibbs- 
Boltzmann equilibrium one. The diffusion constant, when it exists, of the system is 
defined by the late time behaviour of the mean squared displacement as 

(X^) ~ 2DKet. (10) 

Therefore using this notation the bare diffusion constant of the particle in the absence 
of the field ip will always be given hj = kq = 1. 
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We will find here that the low temperature phase is characterised by an anomalous 
sub- diffusive behaviour 

(x?)~^'^ (11) 

where the exponent associated with the anomalous diffusion u < 1/2. We then argue 
that the problem can be mapped onto an effective trap model which allows us to predict 
the value of u. This prediction, in one dimension, is confirmed by a formal replica 
calculation applied to a first passage time problem. It is also supported by our numerical 
results. 



2. Exact Results 



The diffusion constant for a particle in a one dimensional random potential, which is 
statistically invariant under translation, can be obtained exactly in a number of ways 
[T^ ITHj . It is given by 

_ 1 

~ {exp{md{exp{-md ^^^^ 
where the subscript d indicates the disorder average over the field taken at any point. 
An interesting point about this formula is that the resulting Kg is independent of the 
sign of f3. In the case where ip = 0^/2 where the one-point distribution function for 
is Gaussian and given by 

p(0) = ^exp(-0V2) (13) 



/27r 

we have the result 



ne = ^/l-/3' (14) 

and so at /?g = 1 a dynamical transition thus occurs where the diffusion constant 
becomes zero. This transition is really present, and below this transition temperature 
the diffusion will become anomalous. 

Note that if the energy barrier corresponding to a coarse grained region is given by 
^ with the distribution of given by equation (fT^. then the induced distribution on 
the trapping time r given by equation ^ is 

1 

p{r) = . 1 (15) 

r^+^^7ra/31n(^) 

The above trapping time distribution is almost Levy like apart from the logarithmic 
term. Below the critical temperature given by ajSc = 1 we see that the mean value 
diverges, as can be seen directly. The trap model has been extensively studied in finite 
dimensions in this regime. The precise definition is as follows, the system consists of a 
series of traps on a finite dimensional lattice. Associated with each site i is an trap whose 
residence time, the time before the next jump is made, is exponential with quenched 
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mean Tj. If the distribution of the are independent and identically distributed with 
distribution 

P(-) - (16) 

then we see that diverges for /i < 1. For > 1 in one dimension we have that 

(X2)~Ct ;/i>l (17) 

which is normal diffusion. In the anormal phase, < 1 the mean squared displacement 
behaves as [inj 

{Xf) ^Ct^ ]D = l 

~ C t'^(ln(t))i-^ ;D = 2 

~ C ;D > 2. (18) 

From equation (fT3j) we see that up, to logarithmic corrections, the model studied 
here should correspond to a trap model where the exponent of this anomalous diffusion 
is given by 

^i = \ (19) 

The low temperature phase of the model should thus be characterised by anomalous 
diffusion. In fact, in what follows, we shall see, both from analytic arguments and 
numerics, that the effective value of a in all the above is in fact 1. The difficulty of trap 
models with quenched disorder is that a particle may in general visit a trap several times 
and thus it is not necessarily a good approximation to draw a new Xj from the quenched 
distribution each time an already visited site is revisited. This approximation is called 
the annealed approximation and fails badly in one and two dimensions where the random 
walk is recurrent The model corresponding to the annealed approximation is the 
annealed model where each time a particle visits a site the trapping time is redrawn 
from the distribution of waiting times (which is site independent). Above the critical 
value /ic = 1 the annealed model in one dimension diffuses as 

(X?) ~ C t + 0{y/t) ;/i > 2. (20) 

However, still above = 1 but with 1 < < 2 we find 

(X,2) ~ C t + 0(tM) ;1< ^ < 2. (21) 

This change in the exponent of the sub-leading temporal behaviour stems from the 
fact that the variance of the time spent in a trap rf — Ti^ diverges, while the mean 
value remains finite. Thus we see that there are two regimes of normal diffusion, a 
high temperature one corresponding to /i > 2 as given by equation (j2n|l and a low 
temperature normal phase with 1 < /i < 2 and where the corrections to pure diffusion 
change from 0{\/t) to 0{t'^). The key point is the divergence of the variance of the 
trapping time at a given site and we thus expect a similar transition to occur in the 
quenched version of the model. For ji <1 the annealed model gives 

(X,2) ~ Ct^" (22) 
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in all dimensions. Thus we see that, in low dimensions, the exponent associated with 
the anomalous diffusion in the annealed model is not the same as that associated with 
the quenched model. This is precisely due to the recurrence of random walks in two and 
less dimensions. 

For our later comparison with numerical simulations we will need to compute the 
correction to normal diffusion in the high temperature limit. As random walks are 
recurrent in one dimension transient effects are very important and unless one can 
predict them the extraction of a diffusion constant from a numerical simulation becomes 
extremely difficult. 

In the high temperature phase we write 

(X^) ~ 2/tet + 0(t^). (23) 

Now if the corresponding trap model is characterised by a = 1 then we have fi = 1/(3. 
Adapting the scaling argument found in ^H] to the finite time corrections we predict 
that 

e = 3/4 for /i > 2 

= 1 + ^ = (1±1 for 1< /X < 2, (24) 
2 2/i 2 ' V ; 

The above results for the exponent 6 turn out to be in good agreement with our numerical 
results. This supports both our scaling arguments and the identification of our model 
with a trap model characterised by the exponent ^ = 1/(3 (equivalently a = 1). The 
prediction of equation shows that the finite time corrections to diffusion are always 
important, and as one approaches the transition aX (3 = 1 they become of the same 
order as the normal diffusion term. 

Returning to the Langevin problem one may be tempted to argue that the transition 
seen in one dimension could be rather pathological in that the system is obliged to 
overcome all energy barriers in one dimension. Unfortunately no general exact results 
exist in higher dimensions, however in two dimensions if the field i() is statistically the 
same as —%() a duality argument can be used ^7j (in a rather indirect way) to show that 
_ 1 

" (exp(-/5^))/ ^^^^ 
This clearly allows us to compute Ke in two dimensions in the case where %() Gaussian. 
We can however get an exact result for another potential by judiciously choosing a field 
of the form 

where (p and 0' are independent Gaussian fields with the same statistics. Clearly i/j ~ —ip 
in the statistical sense and we thus obtain the exact result 



1 - /5^ (27) 

showing that the same dynamical transition occurs in two dimensions for this particular 
choice of ip. 
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If one wants to simulate the Langevin process considered here on a computer, 
without resorting to a lattice model, one may use a standard method to generate a 
Gaussian random field fH], where one takes 



N 



cos(k„ ■ X + e) 



(28) 



n=l 



where N is the number of modes chosen. When, as we assume here, A(0) = 1 one takes 
each wave vector k„ independently from the distribution 



^(k) = -7^^ J rfx exp(ik ■ x)A(x). 



(27r) 

For instance here we will carry out simulations with 

^2 



A(x) 



exp^-yj 



(29) 



(30) 



and so each component of a vector k„ is a Gaussian of mean zero and variance one. 
In the limit of large N the field so generated becomes Gaussian but it is interesting to 
see what the effect of a finite number of modes is from both a theoretical and practical 
(for comparison with numerical simulations) point of view. We start with the simple 
Gaussian case where we have 



(exp 



1 / r2~ 



cos(e) 



TV 




(31) 



where In is a modified Bessel function llQii. The function In has the small z behaviour 



given by [T 

(1^2)2 

and the large z asymptotic behaviour [1 



+ 0(/ 



exp(2;j 



V2 



TTZ 



l + 0(- 

z' 



(32) 



(33) 



From equation (jH^ we thus find that in the limit N ^ oo but keeping f3 finite we find 



Kg = exp(— 



exp(-— ) 



(34) 



in one dimension. We thus obtain the super- Arrhenius law associated with as expected. 
However at much lower temperatures such that j3 ^ A^^ we use the asymptotic form 
equation to obtain 



N 



exp (-2/?V2iV) ~ exp 



2V2iV. 
7^' 



(35) 



and thus we recover the Arrhenius law at extremely low temperatures. Note that the 
term AE in equation ()35p is simply the maximal difference in energy possible for the 
field expressed as a finite Fourier series as in equation PH|) . 
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As stated earlier when ip = (j) we have the statistical equivalence necessary to use 
the result equation (^3)) and we find the corresponding two dimensional results 



for 



Ke = exp(-y) = 

oo and f3 finite , and 



exp( 



2T2' 



N 



exp 



'2N 



exp( 



T 



(36) 



(37) 



for (3 S> N"^. We see again that the Arrhenius law is recovered in this limit, but in 
contrast with the one dimensional case, the term AE in equation ()37j) is only half of the 
maximal energy difference between two points. This makes sense as in two dimensions 
one can go around this maximal energy barrier. 

The case of a squared Gaussian potential is treated similarly using 



(exp(/3V^))d 



dz 

dz 
sp2^ 



exp( 



exp( 





— r duexpiNfiu)) 



where 



u 



/(«) = --+ In 



h ud2(3 



(38) 



(39) 



The integral in equation (jHSj) can now be treated in the saddle point approximation. 
From the small z expansion of Jq of equation we see that, about m = 0, f{u) takes 
the form 



f{u) 



V 2 iQ V ; 



(40) 



We therefore see that the dynamical transition occurring at /? = 1 is mathematically 
equivalent to a mean field ferromagnetic transition ! Note that if (3 is negative then no 
transition occurs so the term {exp{— jSip)) behaves analytically as A^ ^ cxo. In the high 
temperature phase /3 < 1 therefore we have 



(exp(/3V'))d 



1 



1 



A^(l 



^ + 0(J-^ 



which gives 



8N {1- (3^y ^m' 



(41) 



(42) 



When (3 > 1 the saddle point is no longer at -u = and the function / has a maximum 
value greater than zero at the points ±Uc where the maximum is attained. Here we find 
that 

2 



;v/(l + /5)|/"(wc)|exp {-Nf{u,)) 



(43) 
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For /5 ~ 1 we find that Uc = 2^^/^/ — l/P and consequently 




(44) 



At low temperatures we can use the asymptotic behaviour in equation (jH^ to obtain 
Uc ~ v^2^ and thus 



Again note that that the last part of the above equation indicates an Arrhenius law 
corresponding to the maximum difference in energy between two points, the lowest 
value of ip being zero and the maximum being A^. In two dimensions for the symmetric 
potential of equation (pUj) we find exactly the same results as for the one dimensional 
potential. 

To give the reader a feel for what happens when the number of modes is finite we 
have plotted the N = oo result, equation (fTlj) . against the corresponding results for 
A^ = 64 and A^ = 128, often considered to be sufficient for simulation purposes, which 
are evaluated numerically using the exact equation (jHHj) . We see that even for A^ = 128 
that the value of is still of order 1 ! 

Given the above, it is interesting to examine the case of a diffusing dipole as studied 
in [13]. In the case of positive temperature the potential ip = — |(V0)^ has localised 
minima of varying depths which can be regarded as traps. Therefore in this case one 
would think that the trap model picture would be a reasonably good approximation. 
If we take the equation © with a = 1 and apply it to this case (authors of [TH] used 
A(x) = exp(— x^/2)) one finds that the transition occurs at /? = 1. Interestingly the 
numerical results of ^31 indicate a cross over from a concave behaviour to a convex 
behaviour of Ke in this region, rather reminiscent of the behaviour of the finite A^ values 
in figure near the dynamical transition. It is also interesting to note that the third 
and fourth order perturbation theory results of ^31 give values of Kg which vanish close 
to (3 = 1 (extremely close in the case of the third order result). Having said this we 
have no rigorous proof that there should be a transition in the model studied in |lHj . 
However, for the reasons presented here, a finite size scaling analysis could resolve this 
issue. In ^3] there is no evidence for a transition at negative temperatures. 

To summarise, our exact results have shown that, in the squared Gaussian cases 
studied here, the high temperature regime is characterised by a relaxation time r (using 
the fact that the correlation length scale is 0(1)) which is finite but diverges as T —>■ Tg 
as 



However when the potential is constructed from a finite number of modes we find that 




(45) 



1 



(46) 




for T < Tg. 




(47) 
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Figure 1. Exact value of the diffusion constant in one dimension for a squared 
potential. The solid line shows the exact Gaussian result {N = oo modes) and the 
dotted and dashed lines are for N = 128 and TV = 64 modes respectively. Note that at 
the Gaussian transition temperature given hy (3 — 1 the corresponding Ke for N — 128 
modes is Ke(l) ~ 0-3 



where NAe{T) can be interpreted as a temperature dependent energy barrier. In one 
dimension this barrier height tends to the maximal difference in potential possible 
between two points as the temperature approaches zero and in two dimensions it tends 
to half the maximal difference. 

In one dimension one can obtain information about transport behaviour by 
analysing first passage time. This method can be used to calculate the diffusion constant 
in a one dimensional potential jHj and was used in to determine the exponents of the 
anomalous diffusion in a Gaussian potential with long range (logarithmic) correlations. 
To simplify the calculation one considers diffusion on the half line with reflecting 
boundary conditions at the origin 0. We denote by T{L) the average (in the thermal 
sense of averaging over realisations of the white noise) time at which the tracer particle 
starting at the origin first reaches the point L. This average first passage time is given 
by (see [20] and references within) the following double integral which depends on the 
realisation of the driving field ip 




(48) 
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In general we can compute the moments of T by rephcating the double integrals to yield 
upon disorder averaging 

rL 



T{LY = / n dyadxa9{ya - Xa) det (J + (3AK) ^ , 

where the matrices A and A are symmetric 2n x 2n matrices with 
Aab = A{ya- yb) 

and 

Aab = - ^ab for 1 < a, 6 < n 

= forl<a<n, n + l<6<2n 
= 5ab for n + 1 < a, 6 < 2n, 



(49) 
(50) 



(51) 



where we have used the notation Xa = yn+a- For the first moment n = 1 we find, where 
it is finite, that 



1 



T(L)= Y\dydxe(y-x) 

^0 {1 - f3^ + (3^A^{y - x))^ 



(52) 



For /3 < 1 we may write 



r(L) 



1 



2L2v/r^^ Jo 
which becomes for large L 

T(L) = - LC 

where 



f / dydx 9{y—x) 



(l + ^A2(^- 



2v/r^ 



c 



dz 



1 



1 



,(53) 



(54) 



(55) 



; + ^A2(z))^ 

when it is finite. The leading order in L of the right hand side above yields the correct 
high temperature diffusion constant. 

Inspection of the double integral in equation (fS^ shows that, assuming that A(x) is 
monotonically decreasing, the value of T(L) is finite as long as A(L)2 > (/3^ — 1)//3^. This 
is clearly always the case in the high temperature phase (3 <1. In the low temperature 
phase there is now a length scale L^, which is temperature dependent: for L < 
T{L) is finite, and for L > Lc it is divergent. For a monotone A this means that 
Lc decreases with temperature and we may tentatively relate this with the idea of the 
system progressively freezing on smaller and smaller length scales as the temperature 
is decreased. This image of the system becoming frozen on smaller and smaller length 
scales is one which has been used extensively to interpret experiments on spin glasses 

We will now present an argument which will give a prediction for the exponent of 
anomalous diffusion in the low temperature phase. The argument is similar in spirit to 
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that of j2n| for long range Gaussian potentials in one dimension, although it is difficult 
to make it as rigorous as in the long range case. First notice that as one approaches 
the transition temperature the terms causing the divergence in equation (jH^ are those 
where x and y are far away from each other and hence uncorrelated. Secondly note 
that in the equation PH|) we expect (in the case of > 0) that the first passage time is 
dominated by the maximal value of (p. Let us denote this value by y^. In the replicated 
averaged formula equation (jl^ we thus expect the most divergent term to occur at 
Va = y* for all a. Also we expect that their are many values of x which contribute 
to the most divergent term in equation ()49|) and these are points where 0(x) = 0. In 
general theses points will be well separated and far from y^,. Thus over most of the 
replicated interval \xa — xi,\ ^ 1 and \xa — y*\ ^1. This means that one can make the 
approximation in the integral that A{xa — Xb) ~ and A{xa — y*) ~ and thus the 
replicated integral in equation ()49j) can be approximated as 

T(L)" ^ dy, r dxa det (/ + /?SA)"5 (56) 

JO JO 

1 

« det(/ + /55A)"" 



where 



n + 1 

Bab = — 1 for 1 < a, 6 < n 



= for 1 < a < n, + 1 < 6 < 2n 

= 5ab for n + 1 < a, 6 < 2n, (57) 

The eigenvalues of B are easily calculated, there being 1 equal to —n(3, n — 1 equal to 
and n equal to 1. 
This now gives 

TWr ~ (58) 

(l + n)(l- 71/3)5(1 + /5)t 

Now if we choose the exponent n to be very close to we find that the equation 
(|KHjl is indeed diverging and should be the dominant contribution. This yields 



T(L)^~L^+^ (59) 

and thus dimensionally we have 

(Xf) ~ t^ (60) 

in agreement with the arguments in ^l] and references therein. 

In the presence of an external field the potential becomes ^ip"^ — hx. The same line 
of reasoning now gives for large L 

T{LY dy^ j^" dXa exp (^-nhy^ + ^ hx}j det (/ + pBA)'^ (61) 
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Again we set n very close to ^ 



and dimensionally find the form 



(62) 



again in agreement witli tlie results of [Hj, not only for the temporal exponent but also 
for the exponent associated with h. The above calculation also show that the value of 
a associated with the effective trap model defined by equation ^ for this problem is 
indeed a = 1. 

3. Numerical Simulations 

In this section we will verify the analytical results of the previous sections via two distinct 
types of numerical simulations. We begin with a direct simulation of the Langevin 
process using a second order Runga-Kutta integration of the Langevin equation as 
developed in references j22j and generating the random field via equation p8|l as 
proposed originally in [TH]. This simulation technique allows access to the asymptotic 
regime, where the relevant transport coefficients can be evaluated, at high temperatures 
in the regime of normal diffusion and at temperatures below Tg, but not too low. Clearly 
in order to attain the asymptotic regime at low temperatures we need to diffuse over 
a sufficiently large distance. However as the diffusion becomes progressively slower the 
times needed to attain this regime become prohibitive from the computational point 
of view. We expect that the low temperature phase should be described by a one 
dimensional trap model but, as discussed above, it is not completely obvious what 
parameters we should take for this trap model. We have therefore constructed an 
effective trap model directly from a realisation of the random field as generated by 
equation (^Hj). The procedure used is the following. We take a realisation of the field ip 
over a suitably large interval. Each minima a of the field is associated with a point a 
on the line at its actual position denoted by Xa- We then calculate the energy barrier to 
move to the nearest minima to the left a—1 and the right a+l. These two energy barriers 
are denoted by AE^^ and AE^^ respectively. The particle then is taken to jump to the 
left /right at an exponentially distributed random time T^/ij with average value given by 
the Arrhenius law tl/r = exp (^(3AEj^^-^^^ (the over all scaling of time is unimportant 
to determine the exponent associated with the anomalous diffusion). From the site a 
Tl/r are generated numerically and the particle hops to the site corresponding to the 
shorter time. The time of the simulation is then increased by this shorter time. In this 
way we can assure that the system has diffused sufficiently far and the total simulation 
time is independent of the physical time, allowing us to attain the asymptotic regime. 

Our numerical integration of the stochastic differential in the high temperature 
normal phase yields the following results. Shown in figure is the value of the diffusion 
constant Kg determined from a fit (X^) = 2Kf.t + Bt^ , where 9 is the exponent predicted 
by equation ()24|). The numerical result is generated by averaging over 10000 particles in 
each realisation of the random field and averaging over 100 realisations generated using 
128 modes. The total time of the simulation consisted of 5000 integration steps using 
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Figure 2. Fit of simulation results (circles), assuming an exponent 9 given by equation 
(|24|l . of the diffusion constant in one dimension using N = 128 modes. The dashed 
line shows the analytic result for = 128 modes and the solid line that in the limit 
N ^ oo. 

At = 0.5 in the second order stochastic Runga-Kutta. Again we see that the above fit 
is in good agreement with the theoretical result for = 128 (dashed line) modes. Also 
shown is the A^ — > oo result (solid line) 

In figure (0) we show the fitted value of the exponent 9 assuming the exact result 
(for A^ = 128 modes) for k^. We see that the result is in good accord with the scaling 
prediction of equation for 9. 

For temperatures not too far below Tg we may use a direct numerical integration 
to estimate the exponent u. The results of these simulations are shown in figure (j3]) 
(empty squares). We see that the numerically measured value of the exponent is close 
to the predicted one up to /3 ~ 1.5 but after it departs from the predicted value. We 
believe that this because we are not carrying out the simulation for sufficiently long 
times. Indeed when we simulate the trap model, where simulation time is no longer a 
problem, we see from figure (j3)) (filled circles) that the agreement with the theoretical 
result is much improved for the larger values of 

4. Conclusions 

We have studied the dynamics of a particle diffusing in an potential which is given by the 
square of a Gaussian potential whose statistics are translationally invariant in space and 
whose correlation function is short ranged. In contrast to the Gaussian case, there really 
is a dynamical transition for this model. The transition manifests itself as a crossover 
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Figure 3. Fit of simulation results (circles), of the exponent 6, using the analytic 
result for Kg in one dimension with N = 128 modes. The solid line shows the scaling 
prediction equation for 6. 
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Figure 4. Fit of exponent of anomalous diffusion i/ obtained from the one dimensional 
effective trap model (circles) and direct integration of the stochastic differential 
equation (squares) shown against the prediction v = 1/(1 + /S) (solid line). 
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between a high temperature diffusive regime and a low temperature regime where the 
particle diffusion is anomalous and more specifically is sub-diffusive. We showed how 
the diffusion constant could be computed in one dimension and in a special case in two 
dimensions. Interestingly when the Gaussian field is constructed using a finite Fourier 
series we see that, instead of vanishing at the dynamical transition temperature, the 
diffusion constant obeys an Arrhenius form dependent on the maximal energy barrier 
present in the system. Explicitly we have shown 

Cexp(-^), (63) 

at low temperatures. This result is physically intuitive and our calculation allows for a 
complete understanding of the finite size scaling underlying this dynamical transition. 
We also showed that these finite scaling effects can be important for the number of 
modes typically used in numerical simulations. In addition, in one dimension wc show 
that the low temperature phase can be described in terms of a trap model where the 
energy barriers are assumed to have the same statistics as the energy function itself. 
The resulting sub-diffusive behaviour for the low temperature regime in one dimension 
is 

{X^)^t^. (64) 

This scaling, form and that for the mean displacement in the presence of a uniform drift 
was also obtained via a, non-rigorous, replica based computation of the moments of the 
mean first passage time. Finally in the case of a finite number of modes it is a natural 
question ask which time scale must one go to in order to see the finite size effect diffusive 
behaviour in the low temperature regime as opposed to anomalous diffusion. A simple 
way to estimate this cross over time tc is to equate the mean squared displacements in 
the anomalous and diffusive cases, i.e. 

tf' ~ exp(-^) t,, (65) 

which yields at low temperatures 

tc ~ exp (^^^ . (66) 
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